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ABSTRACT 

We investigate the relationship between galaxies and metal-line absorption systems in a large-scale 
cosmological simulation with galaxy formation. Our detailed treatment of metal enrichment and 
non-equilibrium calculation of oxygen species allow us, for the first time, to carry out quantitative 
calculations of the cross-correlations between galaxies and O VI absorbers. We find the following: 
(1) The cross-correlation strength depends weakly on the absorption strength but strongly on the 
luminosity of the galaxy. (2) The correlation distance increases monotonically with luminosity from 
~ 0.5 — l/i _1 Mpc for 0.1-L, galaxies to ~ 3 — 5/i _1 Mpc for galaxies. (3) The correlation distance 
has a complicated dependence on absorber strength, with a luminosity-dependent peak. (4) Only 
15% of O VI absorbers lie near > L Zj * galaxies. The remaining 85%, then, must arise "near" lower- 
luminosity galaxies, though, the positions of those galaxies is not well-correlated with the absorbers. 
This may point to pollution of intergalactic gas predominantly by smaller galaxies. (5) There is a 
subtle trend that for > 0.5L Zj * galaxies, there is a positive correlation between absorber strength and 
galaxy luminosity in the sense that stronger absorbers have a slightly higher probability of finding 
such a large galaxy at a given projection distance. For less luminous galaxies, there seems to be a 
negative correlation between luminosity and absorber strength. 

Subject headings: stars: abundances — supernovae: general — galaxies: formation — cosmology: 
theory 



1. INTRODUCTION 

Cosmological hydrodynamic simulations have 
shown that most of the s o-call ed missing baryons 
(jFukugita. Hogan. fc Peebles! Il998f) may be in a fila- 
mentary network of Warm-Hot Intergalactic Medium 
(WHIM; ICen fe OstrikerlH99l iDave et aLllMiT ). Visual 
inspection of simulations suggests that the WHIM is 
spatially correlated with galaxies. This is consistent 
with the observed large-scale structure of galaxies, as 
well as the physical expectation that both galaxies and 
intergalactic medium (IGM) are subject to the dominant 
gravitational force of dark m atter which tend s to lead 
to such large-scale structures (jZerdovichlll970h . 

Moreover, the WHIM may provide a primary con- 
duit for matter and energy exchanges between galax- 
ies and the IGM. Thus, a detailed understanding of the 
WHIM may shed useful light on galaxy formation (e.g., 
iCen. Nagamine. fc Ostriken [20051) . The O VI AA1032, 
1038 absorption line doublet in the spectra of low-redshift 
QSOs provides a valuable probe of the WHIM (e.g.. 
Leaner et all 120071 iShull fc Danforthl [20071 iTripp et al l 
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simulations to make predictions of the correlations be- 
tween O VI absorbers and galaxies, which may be used 
for detailed comparisons with upcoming Hubble Space 
Telescope Cosmic Origins Spectrograph observations. 
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We use the simul ation from lCen fc Ostrikerl (|2006l ) and 
I Cen fc Fan g (2006), to make predictions of the relation- 
ship between O VI WHIM absorption and the presence 
of galaxies. This cold dark matter simulation assumes: 
M = 0.31, fl b = 0.048, Q A = 0.69, a$ = 0.89, H = 
100ft- = 69 km s Mpc , n s — 0.97, co-moving box size 
85/i —1 Mpc and 1024 3 cells. The simulation follows star 
formation using a physically-motivated prescription and 
includes feedback processes from star formation to the 
IGM in the form of UV radiation and gal actic superwinds 
carry ing energy and metal-enriched gas (|Cen fc Ostrikerl 
2006). In star- formation sites, "star particles" are pro- 
duced at each time step, which typically have a mass of 
10 6 M Q . Galaxies are identified post facto using the HOP 
grouping scheme (|Eisenstein fc Hutlll998l see Nagamine 
et al. 2001 for details) on these particles. This group- 
ing scheme provides a catalog of galaxies containing 3D 
positions, peculiar velocities, stellar masses and ages. 
The catalog consists of 33,8 87 galaxies. Stellar popu- 
lation synthesis models from iBruzual fc Chariot] (2003) 
are used to compute luminosities in all SDSS bands. For 
the purposes of this paper, we focus on the Sloan z-band, 
which is centered at 9100 A, and has a width of 1200 A 
(jFukugita et al.l[T996h . For galaxies, the flux in this band 
is most tightly correlated ( of all the Sloan bands) with 
the total stellar mass (e.g- lKauffmann et al.l 12003 ). We 
find that the luminosities of field galaxies in this simu- 
lation follow a Schechter function down to ~ 10~ 35 Z/ Zj * 
(L Z: * = 3x 10 11 L Q , Ganguly et al. 2008, in preparation). 

We extract a 100 x 100 grid of sight-lines uniformly 
separated by 850/t" 1 kpc. For each sight-line, a synthetic 
spectrum of O VI A1032 absorption is generated, taking 
into account effects due to peculiar velocities and thermal 
broadening. For each sight-line, we decompose the spec- 
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Fig. 1. — We show the cumulative fraction of O VI absorbers that 
have galaxy at least as luminous as L z /L z>t within a cylindrical 
volume centered on the absorber. We show results for cylinders of 
two different radii, and a velocity depth of 2000 km s — . In the 
top panel, we show different cuts in the integrated O VI column 
density. In the bottom panel, we show cuts in A1032 equivalent 
width. In each panel, a horizontal line is drawn where 50% of the 
absorbers are accounted for. A vertical line is drawn where each of 
the curves crosses this fiducial. Black dashed lines show the results 
for uncorrelated absorber and galaxy positions. 

trum into individual Gaussian-broa dened component s 
with an algorithm similar to AUTOVP (|Dave et al.lll997f) . 
Finally, we associate grouped components into systems; 
this is important since, physically, it does not serve our 
purpose to treat individual components separately when 
comparing to the locations of galaxies. [A single physical 
system like a galactic disk/halo is often composed of mul- 
tiple components. This is a result of complicated velocity 
structures and the dumpiness of gas. Historically, older 
surveys did not have the spectral dispersion needed to 
resolve individual components. Therefore, it is desirable 
to re-group these components back to physically inde- 
pendent systems.] We accomplish this by computing the 
one-dimensional two-point correlation function of com- 
ponents and identify a characteristic velocity scale. We 
find that, on small scales, absorption-line components 
are correlated out to a velocity separation of —±300 km 
s _1 , which we adopt as the characteristic velocity inter- 
val to identify systems. For each system, we record the 
integrated flux- weighted centroid redshift, O VI column 
density, and A1032 rest-frame equivalent width. A total 
of —180,000 components were found to be grouped into 
-21,000 systems with W A (A1032) > 1mA. 



3. ANALYSIS & RESULTS 

We consider the probability that an O VI absorber of a 
given equivalent width or column density lies within some 
distance of a galaxy with a certain luminosity. We take 
three approaches to address this question, but we must 
first tackle the problem of assigning galaxies to O VI ab- 
sorbers. Due to peculiar velocity effects, the exact 3D 
separation between an O VI absorber and a galaxy can- 
not be precisely known, although the projected distance 
in the sky plane can be directly measured. Thus, we 
simply limit absorber-galaxy associations to within the 
physically motivated line-of-sight separation of 1000 km 
s _1 , corresponding approximately to the velocity disper- 
sion of clusters of galaxies. 

Approach #1: For each absorber, we find the most lu- 
minous galaxy within a projected distance of lh _1 Mpc 
or 3h _1 Mpc (and within the aforementioned velocity 
separation). Figure [1] shows the cumulative fraction of 
O VI absorbers that have a galaxy of luminosity > L z 
within this cylindrical volume. We show this distribution 
for O VI absorbers of different integrated column density 
cuts (upper panel) or A1032 equivalent width cuts (lower 
panel). For comparison, we also do the same exercise 
for a sample of galaxies placed randomly in the box, but 
with the same luminosity function and the same total 
number of galaxies. 

Approach #2: We relax the projected distance require- 
ment, and find the closest galaxy with luminosity > /I/ Z) * 
(with four different / = 0.03,0.1,0.3, 1). Figure shows 
the cumulative fraction of O VI absorbers as a function of 
projected distance. As with Figure [TJ we show different 
column density (upper panel) and equivalent width cuts 
(lower panel). In addition, we show the different cuts 
in galaxy luminosity, and we repeat the exercise with a 
sample of randomly placed galaxies. Note that Figures Q] 
and [2] are complementary. A vertical cut in Figure [2] 
at lh _1 Mpc or 3h _1 Mpc, with the cumulative fraction 
of absorbers plotted against luminosity reproduces Fig- 
ure Q] 

Approach #3; For each absorber, we determine the 
number of galaxies that lie within the 3h _1 Mpc radius, 
2000 km s _1 deep cylindrical volume. In Table [TJ we 
list the fraction of these absorbers that have more than 
n galaxies (n = 0, 1, ...) in that volume. We do this for 
two subsamples of absorbers, those with A1032 equiva- 
lent width larger than 30 mA (typical detection limit for 
HST/STIS and FUSE spectra) and 10 mA(expected limit 
for HST/COS observations). In addition to the equiva- 
lent width cuts in the absorber sample, we also make 
luminosity cuts in the galaxy sample. Current catalogs 
are typically able to reach galaxy lum inosities of 0.1L* fo r 
statistically interesting volumes (e.g.. lStocke et al.| [2006). 

Taking Figures Q] and [2] together, we first address the 
question of which galaxies, and on what scales, are corre- 
lated (if any) with O VI WHIM absorption. In Figure [Q 
at luminosities fainter than 0.3 — 0.5L Zj „, the curves lie 
below the equivalent curves for randomly-placed galaxies. 
That is, finding a lower luminosity galaxy near an O VI 
absorber is less probable than if galaxies were uncorre- 
lated with the O VI absorbers for D = 1 — 3/i _1 Mpc. 
This merely implies that the correlation distance, de- 
fined to be the distance within which there is positive 
enhancement of pairs compared to random distributions, 
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Fig. 2. — We show the cumulative fraction of O VI absorbers 
as a function of projected distance to the closest galaxy with lu- 
minosity > L z . The absolute velocity separation is required to be 
|Au| < 1000 km s _1 . We show different cuts in O VI column 
density (top panel) and A1032 equivalent width (bottom panel). 
In both panels, four families of curves are shown corresponding to 
different luminosity cuts: L z /L z ,* > 0.03,0.1,0.3 and 1. Black 
lines show the results for uncorrelated absorber and galaxy po- 
sitions. Note that vertical cuts at -D m i n = l,3h _1 Mpc plotted 
against luminosity reproduce Figure [l] 

is smaller than l/i _1 Mpc for these low luminosity galax- 
ies, consistent with results summarized in Table 2 (see 
below). At higher luminosities, the curves are above 
the random galaxies. This implies that galaxies with 
> 0.3L Zj * are correlated with O VI WHIM absorption 
with the correlation distance larger than 3/i^ 1 Mpc, again 
consistent with results summarized in Table 2. This same 
information may be seen, in a different way, in Figure [21 
as the families of curves for fainter galaxies lie below 
equivalent ones for random galaxies, but the converse 
is true for higher luminosity galaxies. A perhaps some- 
what counter-intuitive result in Figure 1 is that an ab- 
sorber is more likely to find a faint galaxy within the 
D = l/i _1 Mpc cylinder, if the galaxies were randomly 
distributed, and ~ 30% of absorbers do not find any 
galaxy to the faintest limit simulated. This is due to 
the fact that the galaxies themselves are more strongly 
clustered than absorbers among themselves or between 
galaxies and absorbers. Therefore, absorbers will have 
a lower probability of finding neighboring galaxies than 
when the latter are randomly distributed, beyond the 
cross-correlation distance. 

An interesting feature immediately visible from Fig- 
ures Q] and [2] is that the association between galaxies and 



TABLE 1 
O VI Absorbers vs. Galaxies 



n 


/ = 1 


0.3 


0.1 


0.03 


0.01 







W A (A1032) >10 


mA, ' 


^6500 absorbers 





0.15 


0.36 


0.56 


0.74 


0.84 


0.99 


1 


0.04 


0.17 


0.40 


0.62 


0.75 


0.98 


2 


0.01 


0.08 


0.28 


0.52 


0.67 


0.96 


3 


0.00 


0.04 


0.20 


0.43 


0.60 


0.94 


1 


0.00 


0.02 


0.14 


0.36 


0.52 


0.92 


5 


0.00 


0.02 


0.10 


0.30 


0.47 


0.89 


10 


0.00 


0.00 


0.02 


0.12 


0.26 


0.75 




VK A (A1032) >30 


mA, r 


^3000 absorbers 





0.15 


0.35 


0.55 


0.74 


0.84 


0.99 


1 


0.05 


0.17 


0.39 


0.61 


0.75 


0.98 


2 


0.01 


0.08 


0.27 


0.51 


0.67 


0.96 


3 


0.01 


0.05 


0.19 


0.42 


0.59 


0.94 


4 


0.00 


0.02 


0.14 


0.35 


0.52 


0.92 


5 


0.00 


0.02 


0.10 


0.29 


0.46 


0.89 


10 


0.00 


0.00 


0.02 


0.11 


0.25 


0.75 



Note. — Column 1 denotes the number of 
galaxies, n, that lie within a cylindrical vol- 
ume of radius 3h — 1 Mpc and depth 2000 km 
s — 1 centered on an O VI absorber. In the 
header row, / denotes the galaxy luminosity 
in units of L z ,» (i.e., L z = /L Z: *). The num- 
bers in columns 2-7 indicate the fraction of 
> W A (A1032) O VI absorbers that have >n 
galaxies with luminosities > /L Zj ». 



TABLE 2 

O VI Absorber - Galaxy Correlation 
Length 



Lz 1 Lz,* 


W x (0 VI A1032) [mA] 




1 10 30 50 


100 


i 


3.23 4.15 4.20 4.18 


4.92 


0.3 


1.17 1.40 1.84 1.94 


1.27 


0.1 


0.59 0.82 0.74 0.57 




0.03 







Note. — Column 1 denotes the luminosity 
of galaxies in units of L Zj „ =3xlO 11 L0. The 
numbers in columns 2-6 indicate the distance 
in Mpc out to which > W x (0 VI A1032) 
absorbers are correlated with > L z galaxies. 
Galaxies are required to lie within 1000 km 
s~ 1 of the absorber. 

O VI absorbers depends weakly on either the O VI col- 
umn density or A1032 equivalent width. This suggests 
that, while overall galaxies and O VI absorbers are cor- 
related on small scales, the physical properties of O VI 
absorbers (e.g., strength, kinematics, number of compo- 
nents) themselves do not display any tight correspon- 
dence with nearby galaxies. It seems that the strengths 
of O VI absorbers do not provide useful indicators for 
the properties of nearby galaxies. Physically, it suggests 
that O VI absorbers may arise in the vicinities of galaxies 
in a wide variety of ways through complex feedback and 
thermodynamic processes. Complex interactions involv- 
ing gravity- induced shocks, feedback and photoionization 
appear to have erased or smoothed out any potential 
trend with respect to O VI column density or equivalent 
width. This finding is in accord with observations (e.g., 
iProchaska et~aTl [20061 ). 

From Figure [H it appears that the distance out to 
which O VI absorbers are correlated with galaxies is a 
function of both galaxy luminosity and absorber strength 
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(even though the mere presence of a galaxy is not tightly 
correlated with absorber strength as from Figure Q]). 
Comparison of the black curves in Figure [5J showing 
the results of uncorrelated absorber and galaxy positions, 
with the equivalent families of the curves for the simu- 
lated galaxies and absorbers shows that there is typically 
a distance beyond which it is more probable to find an 
uncorrelated galaxy. This distance changes depending 
on the galaxy luminosity and the absorption strength. 
In the extreme case, there is no distance at which the 
strongest absorbers are correlated with the least lumi- 
nous galaxies. In Table [2j we list the correlation dis- 
tances as a function of galaxy luminosity and absorption 
strength. We note two interesting features from the ta- 
ble: (1) For a given equivalent width limit, the corre- 
lation distance is a monotonic function of the limiting 
luminosity. (2) However, for a given limiting luminosity, 
the correlation distance is not a monotonic function of 
equivalent width limit. This may also point toward the 
eclectic nature of the absorbers as mentioned above. 

It is interesting, however, to consider the number of 
galaxies that may be responsible for producing O VI 
absorption. While the intragroup/intracluster medium 
(ICM) is typically too hot (T ~ 10 6 K) for O VI to 
survive in appreciable quantities, the interfaces between 
warm, denser, photoionized clouds of temperature T ~ 
10 4 K and the ICM are potential locations for O VI pro- 
duction. Examples of such interfaces include the bound- 
aries between Milky Wa y high- velocity clouds and the 
hot G alactic corona fe.g.. lFox et al.|[2005tlSembach et al.1 
120031 ). From Table [fl we find that 85% of O VI, regard- 
less of absorption equivalent width, do not lie near L Zi » 
galaxies. Furthermore, the remaining 15% have at most 
3 nearby L z ,*, comparable the Local Group. This is not 
surprising given that L Zj * are not common. However, 
99% of absorbers do lie near galaxies of lower luminosity, 
even if the presence of those galaxies is not correlated 
over what is expected from randomly placed galaxies. 

Figure [2] shows that <20% of O VI absorbers should 
find an L z ^ galaxy within a projected distance of 
5h _1 Mpc. Of course, one would not likely associate 
O VI absorbers with L z „ galaxies at such large projected 
separations, since, for example, one would have already 
found a nearby > 0.03i Zi » galaxy within lh. -1 Mpc, or 

> 0.1L Z! * galaxy at closer distance with comparable 
probability. In any case, it seems unlikely that L z ^ galax- 
ies at such remote distances are responsible for creating 
the O VI absorbers. The rapid rise of probability in Fig- 
ure [H from D m i n = to D m [ n ~ lh _1 Mpc from galaxies 

> 0.03— 0.1L Z * may reflect a ubiquitous physical connec- 
tion between O VI absorbers and these relatively small 
galaxies, perhaps a result of galactic superwinds being 
able to transport metals to a distance of < lh _1 Mpc 
from these galaxies. This, however, does not necessarily 
exclude larger galaxies from being able to do the same. 
The slower rise of probability in Figure [5] from _D m i n = 
to -D m in ~ 5h _1 Mpc for galaxies > 0.3 — 1L Z .* may be 
a result of the intrinsic correlation of large galaxies and 
small galaxies on these scales. These more detailed is- 
sues will be examined subsequently elsewhere. Our re- 
sults appear to be in broad agreement with obse rvations 
( , e.g.. lStockelt^[20^ . [2W7flTri P p et al.ll2006D . 

We take a more detailed look at the dependencies of 



probability on the O VI column density or equivalent 
width. Closer examination of the curves in Figure [T] 
(lower- right corner of both panels) and in Figure (solid 
curves, lower-left corner of lower panel) reveals that the 
probability of finding a ~ L z ^ galaxy is higher for O VI 
absorbers with a higher column density equivalent width. 
In particular, N(0 VI)> 10 14 cm~ 2 or W x > 100 mA 
absorbers deviate noticeably from the weaker absorbers. 
This is a reversal of the trends from other parts of those 
figures. For smaller galaxies, e.g., > 0.3L Zj * galaxies (the 
set of dotted curves in the lower panel of Figure [2]) , the 
trend is considerably weaker although still visible. For 
still smaller galaxies, the trend is reversed, with weaker 
absorbers having a higher probability than stronger ones. 
These results seem to suggest that these very strong O VI 
absorbers tend to be produced in richer, high density 
environments where the probability of finding massive 
galaxies is enhanced (we will examine this physical link 
elsewhere). The fact that weaker O VI absorbers have a 
higher probability of finding a galaxy than stronger ab- 
sorbers, for galaxies less luminous than ~ 0.2 — 0.3L Z * 
(Figure [1]), once again suggests that these relatively 
weaker O VI absorbers (W\ < 50 mA) are probably pro- 
duced by galaxies of ~ 0.1L Zi *, not by more luminous 
galaxies. This is consistent with the trend in Figure [2] 
above. 

4. CONCLUSIONS 

We investigate the relationship between galaxies and 
metal-line absorption systems in a large-scale cos- 
mological simulation with galaxy formation included. 
Our detailed treatment of metal enrichment and non- 
equilibrium calculation of oxygen species allow us, for the 
first time, to carry out quantitative calculations of the 
cross-correlations between galaxies and O VI absorbers. 
We examine the the cross-correlations between O VI ab- 
sorbers and galaxies as a function of projection distance, 
with the line-of-sight velocity separation between an ab- 
sorber and a galaxy constrained to within ±1000km/s. 
Here are some major findings: (1) The cross-correlation 
strength depends only weakly on the strength of the ab- 
sorber but strongly on the luminosity of the galaxy. This 
result suggests that O VI absorbers are produced ubiq- 
uitously and their physical/thermal properties and his- 
tory vary widely, presumably due to the combined ef- 
fects of gravitational shocks, feedback, photo-ionization 
and cooling processes. (2) The correlation length, how- 
ever, does depend on both the galaxy luminosity and on 
the absorber strength from ~ 0.5 — l/i Mpc for 0.1L* 
galaxies to ~ 3 — 5/i _1 Mpc for L* galaxies. While the 
dependence on luminosity is monotonic, the dependence 
on limiting equivalent width appears to peak at some 
luminosity-dependent value and then falls. (3) Only 15% 
of O VI absorbers lie near > L z « galaxies. Thus the 
remaining 85% must be produced by gas ejected from 
fainter galaxies. The positions of lower-luminosity galax- 
ies is not well correlated with absorbers (i.e. in com- 
parison with randomly-placed galaxies) . This may point 
toward pollution of intracluster gas by many galaxies, 
rather than a single high-luminosity galaxy. (4) For 
> 0.5L Zi * galaxies, there is a positive correlation between 
absorber strength and galaxy luminosity (stronger ab- 
sorbers have a slightly higher probability if finding such a 
large galaxy at a given projection distance). The reverse 
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seems true for less luminous galaxies. On average, these 
results indicate that very strong O VI absorbers tend to 
be produced in richer, high density environments where 
it is more likely to find massive galaxies. 

The spatial resolution of our simulation is ~ 80ft. _1 kpc. 
While this is adequate for resolving large galaxies, it be- 
comes marginal for galaxies in halos of total mass less 
than ~ 10 n M Q . Therefore, some of the smaller galax- 
ies of luminosities 0.01 — 0.03L Z * may be significantly 
affected and their abundances underestimated. In addi- 
tion, O VI systems that would have been produced from 
these under-resolved galaxies may be absent. As a re- 
sult, one should treat the cross-correlation strength be- 
tween O VI absorbers and the low-luminosity galaxies as 



a lower bound. But we hope that the preliminary results 
presented here will provide a useful framework for com- 
parison with upcoming imaging campaigns of galaxies in 
the field of quasars and spectroscopic observations with 
the Cosmic Origins Spectrograph. 
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